Semiclassical Time Evolution of the Holes from Luttinger Hamiltonian 
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We study the semi-classical motion of holes by exact numerical solution of the Luttinger model. 
The trajectories obtained for the heavy and light holes agree well with the higher order corrections to 
the abelian and the non-abelian adiabatic theories in Ref. Q [S. Murakami et al., Science 301, 1378 
(2003)], respectively. It is found that the hole trajectories contain rapid oscillations reminiscent of 
the "Zitterbewegung" of relativistic electrons. We also comment on the non-conservation of helicity 
of the light holes. 
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The field of spintronics holds the promise of using the 
spin degree of freedom for building low-power integrated 
information processing and storage devices 0, Q • Spin- 
tronics devices also promises to access the intrinsic quan- 
tum regime of transport, paving the path towards quan- 
tum computing. Recently, it has been predicted theo- 
retically that a dissipationless spin current can be in- 
duced by an external DC electric field in a large class 
of p-doped semi-conductors |]|. The dissipationless spin 
current arises from the spin-orbit coupling in semicon- 
ductors and several other groups have shown that it also 
applies to a broader class of models 0, . 

The theory of Ref.Q is based on the adiabatic solu- 
tion to the Luttinger model, which describes holes near 
the top of the fourfold degenerate valence band. It was 
pointed out that the abelian adiabatic approximation ap- 
plies for the heavy-hole (HH) , while the non-abelian adi- 
abatic approximation is required to obtain the correct 
result for the light-hole (LH). The adiabatic approxima- 
tion is generally based on the separation of the light and 
the heavy hole bands. However, at the top of the va- 
lence band, these two bands intersect each other, and it 
is not clear to which extent the adiabatic approximation 
is valid. In this paper, we solve the semi-classical trajec- 
tory for the Luttinger model exactly by numerical inte- 
gration of the Heisenberg equation of motion. We find 
that the full trajectory of the holes consists of two parts, 
a rapidly oscillating part, reminiscent of the "Zitterbe- 
wegung" of a relativistic electron Q , is super-posed on a 
smooth part, which is accurately described by the adia- 
batic theory. The separation of the rapid and the smooth 
parts of the trajectory is also similar to the cyclotron and 
the guiding center motion of a charged particle in an uni- 
form magnetic field and a spatially varying potential. In 
this sense, the adiabatic approximation in the spin-orbit 
coupled systems is similar to the lowest-Landau-level ap- 
proximation in the quantum Hall effect. 

The Luttinger effective Hamiltonian |]| with an d.c. 



electric field E = i? 2 z can be written as 
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where 71, 72 are the valence-band parameters for semi- 
conductor materials. Luttinger Q pointed out that there 
are 16 linearly independent spin matrices which can be 
chosen as E, S x , S y , S z , S 2 , S 2 , \S X , S y }, {S y ,S z }, 
{S Z ,S X }, {S x ,S y — S z }, {S y ,S z — S x }, {S Z ,S X — S y }, 
S" 3 , S$, S" 3 , S x S y S z + S z S y S x . The full set of dynamic 
variables in the theory consists of three position opera- 
tors x, y and z, three momentum operators k x ,k y and 
k Zl and the 16 spin matrices listed above. The Heisen- 
berg equation of motion for the expectation value of 
any operator A is determined by a differential equation 
d (A) /dt = {ill} -1 ([A, H]). The equations of motion for 
the momentum and the position operators are given by 



d 
dt 



(K) 

(ky) 
(*.) 



(2) 



and 

d_ 

dt 



(y) 



2a 



(k x ) 

(ky) 

(k z ) 

2 (k x Sl) + {k y {S x , S y }) + (k z {S Z: S x }) 

(k X {S X , Sy}) + 2 (kyS 2 y ) + (k Z {Sy, S Z }) 

(k x {S z ,S x }) + (k y {Sy, S z }) + 2 (k z S 2 z ) 
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where a = h 2 (^i + 572/2)/2m, b = —h 2 ^/™-, c = —eE z , 
and { } represents the anticommutative relation. The 
equation of motion for the spin operators can be obtained 
straightforwardly, but they are lengthy and will not be 
given explicitly here. 

Thus the evolution of momentum is simply determined 
by Eq. (J2J) , and can be solved trivially analytically. Next 
we numerically solve the equations for the spin operators, 
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which depends only on the solution of the momentum, 
not on the position. Finally, we decompose the mean 
value of the product of the momentum and the spin into 
the products of their mean values in Eq. J2J), and numer- 
ically solve for the position operators. For convenience, 
we can always choose a coordinate frame which make the 
hole's initial momentum have no y-component. 

Time evolution of the heavy-hole: The initial 
state of the HH with hclicity A = 3/2 can expressed 
as 0(0) = W(k(0))( 1, 0, 0, 0) T , where C/ f (k(0)) = 
exp(— i<f>S z ) exp(— i6S v ) is defined in Ref. Q, and k(0) 
is the initial momentum. The initial mean value of any 
operator A is (A(0)) = (tp(0) \ A\ -0(0)). From this defini- 
tion of the initial state, we obtain (x(0)), (k x (0)), (S x (0)), 
({S x (0),S y (0)}) etc. as the initial conditions for Eqs. J3 
|3j) and the spin equations. 

In Fig. 01 we plot the trajectory of the HH as a func- 
tion of time. We clearly see that besides the acceleration 
in the z direction and the uniform velocity motion along 
the x direction, there is a side-way drift along the y di- 
rection, which is responsible for the dissipationless spin 
current. We can compare the trajectories between our 
numerical solution and the result from Ref. lj. The 
abelian adiabatic equations of Ref. describe the over- 
all trend very well. However, we see that there are rapid 
oscillations on the exact numerical curve. The frequency 
of the oscillations increases and the amplitude decreases 
as the time increases. This "Zitterbewegung" effect can 
be obtained from the higher orders of the adiabatic ap- 
proximation theory J8|]. The oscillation on z(t) can't be 
seen clearly because the figure space is limited, but the 
oscillation on x(t) is really small which is analyzed as 
below. 

In order to study the higher orders of the adi- 
abatic approximation, we transform the Hamiltonian 
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(JTJ. We assume the wavefunction has the form 
of |\I>(x, t)) = exp(— ieE z zt/%) |it(k, t)), then substi- 
tute |^(x, t)) into the Schodinger equation, so that 
we get a new time-dependent Schodinger equation 
iTidt \u(k, t)) = H' (t) |w(k, t)) , where the new time- 
dependent effective Hamiltonian H' (t) — ak(t) 2 + 
b(k(t) ■ S) 2 ,where k(i) is determined by Eq.©. In 
the adiabatic approximation, we assume |u(k, t)) — 
J2C x (t)exp(-iJ*e x (t')dt)Ut(k)\X), where |A) repre- 

A 

sents any eigenstate of S z , so U'(\s.) |A) is the instant 
eigenstate of H' {t). H' (t)W(\s) |A) = e x (t)W(k)\\), 
where e X (t) = h ' ^ e su °stitute |«(k, i)) into the 
time-dependent Schodinger equation, so that we get the 
equation of C\(t) is 



±C(t)+B-C(t) = 0, 
where C(t) = ( C| Ci C_i C_a ) T , and 
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= i J Q Ae(t')dt' is the dynamic phase, and 
ei(t') — eii(t') is the energy difference of 



where a 
Ae(i') = 

HH and LH. If we choose the initial state C(0) = 
( 1, 0, 0, J , the adiabatic approximation assumes 
that w C_3 ±|(t) <C Cs.it) w 1 is always satisfied. 

So only one equation remains, ^Ci(t) = ^9e la . We 

can solve it after the approximation that both Ae and 9 
are slowly varying functions of t. Then the first-order 
correction of trajectory is xW = C%(i)Ci{t) ■ e~ la ■ 

2 2 

(|| U(k)iJ^W(k) ||) + h.c. This method is applica- 
ble to the other three kinds of holes, too. So we get 
the unified formulas of the first-order correction on the 
trajectory of any helicity state, 



iW = 



A(2A 2 - l)eE z sm26 
2fc 2 Ae 



(l-cos(^t)), 



, _ A(2A 2 - l)eE z sin9 f Ae 



y ' = 



k 2 Ae Sin( X i} ' (6) 
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FIG. 1: A heavy-hole (A = 3/2) position three-component vs. 
time. The red solid lines are numerical results, and the blue 
dash lines are from the formulas of Ref. Q]. The initial mo- 
mentum is parallel to x-axis, and the electric field is parallel 
to z-axis. These conditions are same in other figures. 



From these formulas, we can see that the frequency 
uj = Ae/h will increase while the amplitude (fc 2 Ae) _1 
will decrease as time increases, as shown in Fig. ^ We 
can evaluate the quantities of frequency and amplitude 
in Fig. ^ which agree with Eq. ijfJJl very well. The 
oscillation on x(t) is small because sin2f9 w 0. 
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FIG. 2: The mean values of spin three-component of a heavy- 
hole (A = 3/2) vs. time. The red solid lines are numerical 
results, and the blue dash lines are \ki / k. 



FIG. 3: A light-hole (A = 1/2) position three-component vs. 
time. The red solid lines are numerical results, the blue dash 
line is from the abelian adiabatic theory, and the black dot 
lines are from the non-abelian adiabatic theory. 



Now let's study the applicability of Eq. JBJ). We have 

used the approximation that both Ae and 6 are slowly 

varying functions of i, which is equivalent to Aedt <C Ae 

and 9dt <C 9. They imply the same result, eE z At <C hk, 
which means that the approximation is valid when the 
electric field has not brought large changes in momentum. 
If we assume E z = 1 X 10 3 V/m, and k — 4 x 10 8 m _1 , we 
get At <C 2000T. So we have enough periods of oscilla- 
tions in which Eq. © is applicable. 

Fig. indicates Si(t) w Xki(t)/k(t), which implies the 
approximate conservation of HH's helicity. This can be 
seen clearly in Fig. (J5J. The oscillations show that the 
semiclassical spin vector always precesses around the mo- 
mentum direction as the momentum changes in an elec- 
tric field. The oscillation can be calculated with the sim- 
ilar method above. The deep reason for the HH's helicity 
conserving is the matrix element representing transition 
between A = ±3/2 is zero. But the LH's helicity isn't 
conserved as shown in the next section. 

Time evolution of the light-hole: When we choose 
the initial state as ^>(0) = U(k(0))^( 0, 1, 0, ) T , Eqs. 

and the spins' equations describe the evolution of 
a LH with helicity A = 1/2. The trajectory is shown 
in Fig. |31 and the evolution of spin is showed in Fig. 

The anomalous shift in y-direction is not as large 
as predicted from the abelian adiabatic theory of Ref. 
Q, and the helicity is no longer as conserved as that of 
HH. However, both the trajectory and the evolution of 
spin can be explained in the non-abelian adiabatic theory 
0, 0] , which properly takes into account the transition 
between the two LH states. 

If we confine the problem in the light hole's space, Eq. 



Q is reduced to 




= 0. 



(7) 



It describes the evolution of two degenerate states. The 
solution is 



C(t) 



cos(0 t - 0o) sm(9 t - 0„ ) 
-sin(0 t -0 o ) cos(0 f -0 o ) 



(7(0), (8) 



where 9t is the the polar angle at the time t. So we can 
get the anomalous shift in y directions 

y ±h (t) = C^t)U(k)-id ky U^0<.)C(t) 



± 



3cos(0 t - 20 o ) - cos(30 t - 26> ) - 2cos0 o 



4fco sin 9q 



(9) 



and the evolution of spin is (S(i)) 
Ct(i)E/(k)SC/t(k)C(£), 

S*M (*) = T[7 sin(0 4 - 20 o ) + \ sin(30 t - 20 o ) 

s vM (t) = o, 



(10) 



v 

S zM (t) =±£ cos(0 t - 29 Q ) - - cos(30 t - 29 )}. 

The results from Eq. JjJ} and (|10|) has been plotted in 
Fig. |3 and Fig. they describe the trends of numer- 
ical curves very well except for the rapid oscillation on 
the numerical curves, which have been explained in the 
previous sections due to higher order corrections to the 
adiabatic theory. 

At last, we obtain the anomalous velocity in y- 
direction, 



3eE z 
Ahk 2 



[wx(fi t - 20 o ) - sin(30 t - 29 )]. (11) 
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FIG. 4: The mean values of spin three-component of a light- 
hole (A = 1/2) vs. time. The red solid lines are numerical 
results, and the black dot lines are from the non-abelian adi- 
abatic theory. 



When t = 0, v y ± i(ti) = A(2A 2 - l)eE z k x0 /(hk$), which 
is just the Eq. (7) of Ref. [y] (where fy is given by Eq. 
(S5) of SOM). Eq. I|llfl represent the anomalous velocity 
at any time. 

Unlike the HH, the LH does not always stay as an 
eigenstate, it will evolute according to Eq. 0. Fig. © 
compares the spins' evolution of HH and LH. Obviously, 
LH's helicity is not as conserved as HH, so LH's spin 
can't be always parallel to its momentum like HH. The 
non-abelian adiabatic theory of Ref. properly takes 
this effect into account. 

The adiabatic condition: Ref. jij raised a criti- 
cism by asking why the anomalous shift in Ref. is 
independent of A2. Actually, if A2 = 0, the anomalous 
shift vanishes because the Hamiltonian degenerates to 
an ordinary one without the spin-orbit coupling, and the 
adiabatic approximation is no longer valid. Below can 
we see explicitly that the adiabatic approximation fails 
when A2 is less than a certain quantity. The condition of 
adiabatic approximation || is 



(H^\i\L,(3) 



V3d8 
2 dt 



\fimeE z sin 9 
4fc 3 ?i 2 7 2 



< 1. (12) 



The condition is better satisfied if E z is smaller and 72 
is larger. The small E z ensures that the time-dependent 
Hamiltonian changes slowly, and the large 72 ensures that 
the energy difference between the HH and LH bands is 
large (Ae = — — 72 ), so the transition probability be- 
tween HH and LH is small. 

In most semiconductors, Eq. I|12|) can be satisfied. For 
example as GaAs, 72 = 1.01, kp ps 8 x 10 8 m _1 , if we 
assume E z — 10 3 V/m, 6q — 90°, we get the condition is 
k 3> 0.02/cf- So only a little part in the middle of Fermi 
Ball doesn't meet the conditions. We can neglect them 
when we integrate the whole Fermi ball. 



Conclusion We have studied the motions of the 
heavy-hole and light-hole in a large class of hole-doped 
semiconductor based on the Luttinger Hamiltonian. The 
trajectory of HH has rapid, small amplitude oscillations, 
which can be explained as the first-order correction on 
the trend described in Ref. Q. The trajectory of LH is 
more complicated and the helicity of the LH is not con- 
served. The non-conservation of the helicity invalidates 
the abelian adiabatic approximation. However, the mo- 
tion of LH can be well explained by the non-abelian adi- 
abatic theory. The excellent agreement between the ex- 
act numerical solution of the Heisenberg equation of mo- 
tion and the adiabatic approximation validates the key 
assumptions leading to the dissipationless spin current, 
and addresses the naive criticism raised in Ref. @. In 
the future, we plan to apply the formalisms developed 
in this paper to study the Luttinger Hamiltonian under 
more general external conditions. 

We thank Andrei Bernevig, D. Culcer and Q. Niu for 
helpful discussion. This work was supported by the NSF 
of China under grant number 10174095 and 90103024, 
the US NSF under grant number DMR-0342832, and the 
US Department of Energy, Office of Basic Energy Sci- 
ences under contract DE-AC03-76SF00515. 



[1] 



[2] 
[3] 
[4] 

[5] 

[6] 

[7] 
[8] 
[9] 



S. Murakami, N. Nagaosa, and S. C. Zhang, Science 301, 

1378 (2003); Supporting Online Material (SOM) text for 

previous paper; S. Murakami, N. Nagaosa, and S. C. 

Zhang, Phys. Rev. B 69, 235206 (2004). 

G. A. Prinz, Science 282, 1660 (1998). 

S. A. Wolf et al, Science 294, 1488 (2001). 

J. Sinova, D. Culcer, Q. Niu, N. A. Sinitsyn, T. Jungwirth, 

and A. H. MacDonald, Phys. Rev. Lett. 92, 126603 (2004). 

S. Q. Shen, M. Ma, X. C. Xie, and F. C. Zhang, Phys. 

Rev. Lett. 92, 256603 (2004). 

E. Schrodinger, Sitzungsb. Preuss. Akad. Wiss. 24, 418 
(1930). 

J. M. Luttinger, Phys. Rev. 102, 1030 (1956). 

C. P. Sun, Phys. Rev. D 41, 1 349 (1990). 

X. D. Wang, and X. G. Zhang, [concLmat/0407699 (2004). 



